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The formation of regular patterns is a well-known phenomenon in condensed matter physics. Sys- 
tems that exhibit pattern formation are typically driven and dissipative with pattern formation 
occurring in the weakly non-linear regime and sometimes even in more strongly non-linear regions 
of parameter space. In the early universe, parametric resonance can drive explosive particle pro- 
duction called preheating. The fields that are populated then decay quantum mechanically if their 
particles are unstable. Thus, during preheating, a driven-dissipative system exists. In this paper, 
we show that a self-coupled inflaton oscillating in its potential at the end of inflation can exhibit 
pattern formation. 
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Much recent work has been done on the topic of pre- 
heating in inflationary cosmology. Preheating is a stage 
of explosive particle production which results from the 
resonant driving of particle modes by an inflaton oscil- 
lating in its potential at the end of inflation . 

In regions of parameter space where parametric res- 
onance is effective, much of the energy of the inflaton 
is transferred to bands of resonant wave modes. This 
energy transfer is non-thermal and can lead to interest- 
ing non-equilibrium behavior. Two examples of the non- 
equilibrium effects that can be produced are non-thermal 
phase transitions Q-Q] and baryogenesis ||[|. The non- 
thermal phase transitions induced during preheating can 
also possibly lead to topological defect formation [^o| |l2| , 
even at energies above the eventual final thermal tem- 
perature. Furthermore, non-linear evolution of the field 
when quantum decay of the resonantly produced parti- 
cles is small leads to a turbulent power-law spectrum of 
density fluctuations 13 19]. 



In this letter, we present a new phenomenon that can 
arise from preheating: pattern formation. It has long 
been known that many condensed matter systems ex- 
hibit pattern format ion^. Examples of pattern forming 
systems which have been studied are ripples on sand 
dunes, cloud streets and a variety of other convective 
systems, chemical reaction-diffusion systems, stellar at- 
mospheres and vibrated granular materials. All of these 
physical systems have two features in common. They are 
all driven in some manner, i.e. energy is input to the sys- 
tem, and they are all dissipative, usually being governed 



by diffusive equations of motion. Typically, patterns are 
formed in these systems in the weakly non-linear regime 
before the energy introduced into the system overwhelms 
the dissipative mechanism. Sometimes, patterns persist 
beyond the weakly non-linear regime as well. 

At the end of inflation, the inflaton <j> is homogeneous 
with small perturbations 5<p imprinted on it due to quan- 
tum fluctuations. In chaotic inflationary models, the in- 
flaton then oscillates about the minimum of its potential, 
giving an effectively time dependent mass to fields with 
which it is coupled. The time dependent mass drives ex- 
ponential growth in the population of bands of particle 
wave modes. Many of the fields into which the infla- 
ton can decay resonantly are also unstable to quantum 
decay. For these reasons, at the end of inflation, we are 
considering fields which are driven, due to resonant parti- 
cle creation, and also dissipative, due to quantum decay. 
Therefore, it makes sense to see if there is a region of 
parameter space where the system is in the weakly non- 
linear regime and exhibits pattern formation. 

In our investigation we consider A0 4 theory with the 
addition of a phenomenological decay term to mimic the 
inflaton's quantum decay. This model without the de- 
cay term has been studied extensively in the literature 
[p|,pp^ p^Jl9[ | , and a similar model including the decay 
term has also been studied |^0| . 

Our field equation is 



<t> + 7^ - V 2 + \cf> 3 







(1) 



*For an extensive review of pattern formation in condensed 
matter systems, see |T^ , |l^| 



where 7 is a decay constant and A is the self-coupling 
of the field. Here, we neglect the expansion of the uni- 
verse, although we will comment on the effect of expan- 
sion later. For our calculations, we rescale: t — > t/^/\cf>o, 
x — ► x/^/XipQ and 4> — > 4>4>o, where c/)q is the value of 
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the inflaton at the end of inflation. This gives us a new 
equation 

4> + T<j> - V 2 <j) + (j) 3 = (2) 

where T = 7/ VA0o • 

It should be noted that pattern formation in the infla- 
ton system is conceptually distinct from condensed mat- 
ter systems for at least two reasons. First, the equations 
we study are wave equations with damping, not diffu- 
sive equations. Secondly, we expect wave patterns to be 
formed while the homogeneous mode decays, therefore 
pattern formation will be a temporary phenomenon, at 
least in the model above in which gravity is neglected. 
This should be considered in contrast to the typical con- 
densed matter system, in which energy is introduced via 
boundary conditions (in a convective system) or by a 
vibrating bed (in a granular material system), and the 
energy input is essentially constant. 

For r = 0, the resonant modes lie in the interval 

I < k 2 < V3. (3) 
For r ^ 0, we can introduce ip = §e^ 1 giving 

^5-(V 2 + ^)^ + e- r V 3 = 0. (4) 

Therefore, for small T, we expect the resonance bands 
to be slightly shifted, because the V 2 term is shifted by 
r 2 /4, and we expect the resonance to diminish slowly 
over time, due to the exponential damping of the po- 
tential with time. It should be noted however that the 
resonance structure of the equation can be quite sensitive 
to changes in the potential p9[ . 

As the effect that we are trying to isolate is non-linear, 
we resort to numerical simulation of the field equation. 
We use a leapfrog code which is second-order accurate 
in time and we use fourth-order spatial differences. We 
simulate the field in two dimensions in a box with peri- 
odic boundary conditions which has 256 grid points per 
dimension. 

For initial conditions, we give small amplitudes (~ 
10~ 3 ) to all resonant modes in the box. It is unneces- 
sary to populate non-resonant modes because they will 
decay while the resonant modes grow. In the initial con- 
dition setup, it is crucial that there be many resonant 
modes in the box. 

Early on, we found spurious pattern formation when 
we set the box size such that the resonant mode wave 
number in the box was 3. For such small wave numbers, 
only two resonant modes exist in the box, one along the x- 
axis, and one along the y-axis, giving rise to a misleading 
square wave pattern. Setting the box size such that the 
resonant wave number in the box is 16 is enough to give 
many different resonant modes in the box, but still have 



good resolution of the wave, so this is the box size we 
used. 

We set out to identify the weakly non-linear regime. 
In |l3| it is shown that the self-coupled inflaton system's 
non-linear time evolution proceeds as follows: First, the 
resonant band amplitude grows. Next, when the ampli- 
tude in the resonant band is high enough for non-linear 
effects to become important, period doubling occurs and 
subsidiary peaks develop in the power spectrum. Further 
peaks then develop and the spectrum broadens and ap- 
proaches a power-law spectrum. We tuned T such that 
the amplitude of the resonant band grew, but little period 
doubling occurred. In this regime, only resonant mode 
wavelengths would exist in the box and they would in- 
teract with each other non-linearly. 

In figure [l] we plot a superposition of the power spec- 
trum at various times during our simulation. It is possi- 
ble to see that for the value T = 0.0035, the system stays 
in the weakly non-linear regime for the entire simulation. 
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FIG. 1. A superposition of power spectra taken from a sim- 
ulation with r = 0.0035. We plot the amplitude of the power 
spectra vs. wave number. Note that the period doubling 
modes are only weakly populated, indicating that the simu- 
lation is in the weakly non-linear regime. 

When wave patterns form, the specific pattern which 
arises is due to the non-linear interaction of the wave 
modes. The amplitude of wave modes separated by dif- 
ferent angles grows at different rates. Modes separated 
by angles with the fastest growing amplitudes dominate 
the solution and form the wave pattern. In plots |2| ||, ^ 
and |5|, we show the evolution of the Fourier transform of 
the inflaton. In the final plot, it is clear that the domi- 
nant modes have picked out a preferred configuration in 
Fourier space. In figure |6| we also plot (f)(x, y) in configu- 
ration space at the time when the pattern has formed. It 
should be commented that the preferred angles in Fourier 
space of the pattern would pick out a dodecahedron, but 
this is not a close-packable structure, therefore the pat- 
tern looks more complicated in configuration space, than 
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in Fourier space. After this time, the field amplitude be- 
gins to decay; the wave pattern remains imprinted, but 
decreases in amplitude. 




FIG. 2. Initial conditions for the simulation (F = 0.0035). 
4>(k x , k y ), the Fourier transform of the inflaton is plotted. The 
central peak is the zero mode, and the surrounding ring is the 
resonant mode populated with small amplitudes, t = 0. Only 
the region of interest is plotted. All modes outside of the 
plotted region have zero amplitude. 




FIG. 3. After resonance begins to boost the amplitude of 
the resonant mode. Notice the brightening (increasing ampli- 
tude) of the resonant modes, t = 60. 
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FIG. 4. Non-linear effects have become important by this 
time causing fragmenting of the resonant mode. Modes with 
dominant growth are taking over, t = 210. 




FIG. 5. The final wave pattern, t = 300. 



The important result here is not that we have found 
the preferred wave pattern of the system. With only a 
discrete set of wavemodes in the box, we cannot be sure 
that some other pattern would not dominate if we were 
able to investigate a continuum of modes. However, it is 
clear that some pattern will be picked out by the system 
in this parameter region. Therefore, we can claim to have 
identified pattern formation as a phenomenon exhibited 
by the damped, self-coupled inflaton. 

It would also seem that pattern formation should exist 
in realizations of parametric resonance and preheating 
other than the system we have analyzed, simply because 
of their driven-dissipative nature. 
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In an expanding universe, after a similar variable and 
field transformation as before, the field equations (in con- 
formal time) become 

(p + aTip - V 2 ip - (dr+ -) w+ tp 3 = 0, (5) 

a 

where here ip = a<f>. Note that there is no Hubble damp- 
ing term because the theory without damping due to 
quantum decay is conformally invariant. (We are study- 
ing other models though, to see whether Hubble damping 
can provide the necessary damping for pattern forma- 
tion.) For the expanding case, we expect pattern forma- 
tion for smaller values of T as the effective dissipation 
grows with the scale factor a. Oscillation of the homoge- 
neous mode should be dominated by the ip 3 term at early 
times before ip has had time to be significantly damped. 
We are currently investigating this case. 




FIG. 6. The final wave pattern in configuration space. 
t = 300. The corners are darkened by small amplitude long 
wavelength modes which modulate the wave pattern. 

For wave patterns to have an impact cosmologically, 
they must persist after they are created. In our toy 
model, the resonance band corresponds to wavelengths 
roughly the size of the Hubble radius at the start of pre- 
heating, about 10~ 25 cm. The size of these fluctuations 
today would then only be about 10 m. However we have 
seen that the patterns are modulated by long wavelength 
modes, and it is this non-linear transfer of power to small 
/c-modes that may be important. Thus the size of the 
Hubble radius at the end of preheating may actually set 
the scale for wave patterns. In any case, as the wave 
patterns are regions of energy density and are present es- 
sentially from the beginning of the universe, they may be 
sites of significant gravitational accretion. As dark mat- 
ter seeds it may be possible for wave patterns to change 
the thermal history of baryonic matter. It is this di- 



rection we intend to take in determining the role wave 
patterns may play in the history of the universe. 
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